Finite dipolar hexagonal columns on piled layers of triangular lattice 
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We have investigated, by the Monte Carlo simulation, spin systems which represent moments of 
arrayed magnetic nanoparticles interacting with each other only by the dipole-dipole interaction. 
, In the present paper we aim the understanding of finite size effects on the magnetic nanoparticles 

arrayed in hexagonal columns cut out from the close-packing structures or from those with uniaxial 
compression. In columns with the genuine close-packing structures, we observe a single vortex state 
which is also observed previously in finite 2-dimensional systems. On the other hand in the system 
i <~| with the inter-layer distance set 1/V2 times of the close-packing one, we found ground states which 

O ' depend on the number of layers. The dependence is induced by a finite size effect and is related to 

a orientation transition in the corresponding bulk system. 



In recent years, systems consisting of arrayed single-domain ferromagnetic nanoparticles receive much attention 
of many researchers because of its possibility as a high storage density In such array systems the dipole-dipole 
interaction plays an important role on their magnetic properties. In particular the dipole-dipole interaction originates 
strong finite size effects. For example the peculiar "from-edge-to-interior freezing" of a finite dipolar cube is found 
in our recent workQ- In the present paper, we aim to clarify finite size effects due to the dipole-dipole interaction 
on magnetic nanoparticles arrayed in hexagonal columns cut out from piled triangular lattice with the close-packing 
'"^J ] structures (FCC and HCP) or with those uniaxially-compressed in a direction of the 3-fold axis. 

We perform simulated annealing of the dipolar columns by the heat-bath Monte Carlo method. We employ the 
O \ well-known Hamiltonian of a dipole-dipole interacting system, 

where Sj's denote classic Heisenbcrg spins which represent magnetic moments of nanoparticles. D is a coupling 
constant and a a lattice constant. Other effects in nanoparticle systems are omitted because we focus only on those 
\Q ' by the dipole-dipole interaction. The number of piled layers in the columns examined, JVl, is from 1 to 15. Each layer 
, contains about 100 sites. 

As well known spin configurations consist of in-plane vortices at low temperatures in finite dipolar systems 

on the triangular lattice. In the present simulation, spins in the ground state of a single-layer system also form an 
in-plane single vortex shown in Fig. ^a). We can explain the origin of the vortex state as below. In the bulk system 
the ground state is in-plane ferromagnetic onc||. In the state, Hamiltonian Q possesses the global symmetry of 
0(2), which corresponds to the global in-plane spin-rotation. Thus the direction of the ferromagnetic order is not 
fixed in the bulk system. In a finite system, on the other hand, the symmetry around an edge is reduced into Zi 
£j . by the anisotropy of the dipole-dipole interaction. On the edge the ferromagnetic order is in the direction parallel 
to the edge. A fixed ferromagnetic domain extends from each edge because of this symmetry reduction. In fact, the 
from-edge-to-interior freezing is observed in the single-layer system. The ferromagnetic domains continuously merge 
with neighboring domains and form the single vortex state as a whole. This symmetry reduction is similar to the one 
5-j ' found in the dipolar cube0. We can also expect it around domain walls observed in large systems^. 

The nature of the single layer is kept in multi-layer columns with the genuine close-packing structures. In fact a 
single vortex state is also observed as a whole as shown in Fig. ^b) . We next examine columns with the inter-layer 
distance set 1/yz times of the close-packing one, i.e., columns on close-packing lattices compressed in the direction 
of their 3-fold axis. In the system of 2 layers, the single vortex state is also observed. However we have found a 
multi-fcrromagnetic-domain structure, which is shown in Fig. [2^a) , when TVl is more than 2. In its ferromagnetic 
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FIG. 1: The schematic image of the single vortex state in the single-layer (a) and multi-layer (b) systems. The solid and dashed 
lines denote the directions of the spins and the domain boundary between ferromagnetic domains, respectively. 



domain, spins align in the direction perpendicular to layers and in the direction anti-parallel to those of neighboring 
domains except for those around top and bottom layers. On the top and bottom layers, spins align so as to connect 
smoothly between those in neighboring domains. Spin configuration as a whole form a convective flow pattern of 
ferromagnetic spin chains between the top and bottom layers as shown in Fig. Efb). In the multi-layer hexagonal 
column, the number of the domain decreases from 6 to 3 as Ni, increases from 3 to 15. 

To interpret the results observed numerically, let us first consider the corresponding bulk system. Its ground state 
is considered to consist of multi- ferromagnetic domain |^. We focus attention on a domain in the ground state. Now 
we assume that the domain is enough large that contributions of boundary energy and of inter-domain interaction 
are negligibly small as compared to that of the intra-domain interaction energy, Ed . Because of the symmetry of 
Hamiltonian Ed is described by magnetization of the domain, 

E U = -J{M 2 X + M 2 ) - AM Z 2 , (2) 

where M x , M y and M z denote the magnetization of the domain with J2 a M 2 = M 2 , J and A anisotropy constants. 
The z-direction is set in the direction of the 3-fold axis. This formula holds in all lattices examined in the present 
paper, i.e., FCC, HCP and those with the compression in the z-direction. In these systems if J > A, the in-plane 
0(2)-symmetric state with M z = appears. If J < A, the out-of-plane state with M z = ±M appears. The bulk 
system is, thereby, expected to exhibit a first order transition when the ratio J/A is change. As known previously, 
J = A is for genuine FCC0. We assume from our observation described above that J ~ A is for genuine HCP and 
that J < A for those compressed 1/V2 times at least. 

Next let us introduce upper and lower boundaries into the bulk single domain. If JVl is 2, i.e., the system consists 
of only two boundary triangular lattice layers, the situation of J > A is expected to hold due to the strong anisotropy 
on the boundaries and the in-plane 0(2)-symmetric state is preferred. When we next introduce the lateral boundary 
with the hexagonal shape, the in-plane single vortex state appears as in the single layer described above. As the 
number of layers increases, the effective fields that spins inside a column feel come closer to those of the bulk single 
domain. In columns consisting of multi-layer on the genuine close-packing lattice, the situation of J > A is kept, 
and the single vortex state as a whole is realized (Fig. ^b)), because J ~ A at most in the bulk single domain. In 
multi-layer systems with the compression, on the other hand, the situation of J < A is expected as we have assumed 
above, and the out-of-plane state is preferred. The consequence is the convective flow pattern state shown in Fig. 
which associates with zero total magnetization. The latter result is a very common property of dipolar systems, 
because of which we have to consider a multi-domain state even in bulk (infinite) systems. 

To conclude, the finite size effect, combined with the peculiar anisotropic nature of the dipolc-dipolc interaction 
induces the orientation change in the ground state of dipolar columns consisting of triangular lattice layers as the 
intcrlayer distance or the number of layers is changed. 
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